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- - 1 We propose the most general modified first-order Hofava-Lifshitz gravity, whose action does not 

' contain time derivatives higher than the second order. The Hamiltonian structure of this theory is 

, studied in all the details in the case of the spatially-fiat FRW space-time, demonstrating many of 

' the features of the general theory. It is shown that, with some plausible assumptions, including the 

, projectability of the lapse function, this model is consistent. As a large class of such theories, the 

'~{ • modified Hoi^ava-Lifshitz F[R) gravity is introduced. The study of its ultraviolet properties shows 

' that its 2: = 3 version seems to be renormalizable in the same way as the original Hoi'ava-Lifshitz 

1/^ . proposal. The Hamiltonian analysis of the modified Hoi^ava-Lifshitz F{R) gravity shows that it is in 

' general a consistent theory. The F{R) gravity action is also studied in the fixed-gauge form, where 

, the appearance of a scalar field is particularly illustrative. Then the spatially-fiat FRW cosmology 

(-H ' for this F{R) gravity is investigated. It is shown that a special choice of parameters for this theory 

, leads to the same equations of motion as in the case of traditional F{R) gravity. Nevertheless, the 

^ • cosmological structure of the modified Hofava-Lifshitz F{R) gravity turns out to be much richer 

U , than for its traditional counterpart. The emergence of multiple de Sitter solutions indicates to the 

i-C ■ possibility of unification of early-time infiation with late-time acceleration within the same model. 

Power-law F(R) theories are also investigated in detail. It is analytically shown that they have 

. a quite rich cosmological structure: early /late-time cosmic acceleration of quintessence, as well 

^ ' as of phantom types. Also it is demonstrated that all the four known types of finite-time future 
singularities may occur in the power-law Hofava-Lifshitz F{R) gravity. Finally, a covariant proposal 

' for (renormalizable) F{R) gravity within the Hofava-Lifshitz spirit is presented. 

^ , PACS numbers: ll.lO.Ef, 95.36.-l-x, 98.80.Cq, 04.50.Kd, ll.25.-w 

cn ■ 
O ; 

O . I. INTRODUCTION 

i> ! . 

. !^ ' Recently, it has become clear that our universe has not only undergone the period of early-time accelerated expansion 
(inflation), but also is currently in the so-called late-time accelerating epoch (dark energy era). An extremely powerful 

J-j , way to describe the early-time inflation and the latc-time acceleration in a unified manner is modified gravity. This 



approach does not require the introduction of new dark components like inflaton and dark energy. The unified 
description of infiation and dark energy is achieved by modifying the gravitational action at the very early universe as 
well as at the very late times (for a review of such models, see [ij). A number of viable modified gravity theories has 
been suggested. Despite some indications to possible connection with string/M-theory such theories remain to be 
mainly phenomenological. It is a challenge to investigate their origin from some (not yet constructed) fundamental 
quantum gravity theory. 

Among the recent attempts to construct a consistent theory of quantum gravity much attention has been paid to 
the quite remarkable Hofava-Lifshitz quantum gravity 3| , which appears to be power-counting renormalizable in four 
dimensions. In this theory the local Lorentz invariance is abandoned, but it is restored as an approximate symmetry 
at low energies. Despite its partial success as a candidate for fundamental theory of gravity, there are a number of 
unresolved problems related to the detailed balance and projectability conditions, consistency, its general relativity 



* Also at Tomsk State Pedagogical University 



2 



(GR) limit, realistic cosmological applications, the relation to other modified gravities, etc. Due to the fact that its 
spatially-flat FRW cosmology is almost the same as in GR, it is difficult to obtain a unified description of the 
early-time inflation with the late-time acceleration in the standard Hof ava-Lifshitz gravity. 

Recently the modified Hof ava-Lifshitz F{R) gravity has been proposed [1]. Such a modification may be easily related 
with the traditional modified gravity approach, but turns out to be much richer in terms of the possible cosmological 
solutions. For instance, the unification of inflation with dark energy seems to be possible in such Hofava-Lifshitz 
gravity due to the presence of multiple de Sitter solutions. Moreover, on the one hand, there is the hope that the 
generalization of Hofava-Lifshitz gravity may lead to new classes of renormalizable quantum gravity. On the other 
hand, one may hope to formulate the dynamical scenario for the Lorentz symmetry violation/restoration, caused by 
the expansion of the universe, in terms of such generalized theory. 

In the present work (section|lI| we propose the most general modifled flrst-order Hof ava-Lifshitz-like theory, without 
higher derivative terms which are normally responsible for the presence of ghosts. The general form of the action in 
the spatially-flat FRW space-time is found, and the Hamiltonian structure of the action is analyzed in section IHII 

As a specific example of such a first-order action we introduce the modified Hofava-Lifshitz F{R) theory which is 
more general than the model of ref. [sj. Nevertheless, its spatially-flat FRW cosmology turns out to be the same as 
for the model (this is not the case for BH solutions, etc). Therefore it also coincides with the conventional F{R) 
spatially-flat cosmology for a speciflc choice of the parameters. The ultraviolet structure of the new Hofava-Lifshitz 
F(R) gravity is carefully investigated. It is shown that such models can have very nice ultraviolet behaviour at z = 2. 
Moreover, for z = 3 a big class of renormalizable models is suggested (section |n|. The Hamiltonian analysis of the 
modified Hofava-Lifshitz F{R) gravity is presented in section IIVI The fixed gauge modified Hofava-Lifshitz F{R) 
gravity is analyzed in section [V] 

Section rvTl is devoted to the investigation of spatially-flat FRW cosmology for power-law F{R) gravity. The general 
equation for the de Sitter solutions is obtained. It acquires an extremely simple form for a special choice of parameters, 
when de Sitter solutions are roots of the equation F = 0. The existence of multiple de Sitter solutions indicates the 
principal possibility of attaining the uniflcation of the early-time inflation with the late-time acceleration in the 
modified Hofava-Lifshitz F{R) gravity. The reconstruction technique is developed for the study of analytical and 
accelerating FRW cosmologies in power-law models. A number of explicit analytical solutions are presented. It is 
shown by explicit examples that some of the quintessence/phantom- like cosmologies may develop the future finite-time 
singularity of all the known four types, precisely in the same way as for traditional dark energy models. The possible 
curing of such singularities could be achieved in a similar way as in the case of traditional modified gravity. Some 
remarks about small corrections to the Newton law are made in section fVI 1 1 A summary and outlook are given in the 
last section rVIIIl In the appendix |X] we propose a covariant F(R) gravity that is quite similar to the corresponding 
Hofava-Lifshitz version but remains to be a covariant theory. It seems that it could also be made renormalizable. 



II. GENERAL ACTION FOR HORAVA-LIFSHITZ-LIKE GRAVITY AND RENORMALIZABILITY 

In this section we propose the essentially most general Hofava-Lifshitz-like gravity action, which does not contain 
derivatives with respect to the time coordinate higher than the second order. Its ultraviolet properties are discussed. 

By using the ADM decomposition 6] (for reviews and mathematical background, see [3]), one can write the metric 
of space-time in the following form: 

ds^ = -N^dt^ + glf (dx' + N'dt) (dx^ + N^dt) , i, j = 1, 2, 3 . (1) 

Here N is called the lapse variable and A* is the shift 3-vector. Then the scalar curvature R has the following form: 

R = K'^K.j ~K'^ + + 2V^ (n^V^n"^ - n^W^n") . (2) 

Here R^^^ is the three-dimensional scalar curvature defined by the metric g^j"^ and Kij is the extrinsic curvature defined 
by 



^ (slf - Vf ^iV, - Vf A^,) , K = K\. (3) 

n'^ is the unit vector perpendicular to the three-dimensional space-like hypersurface Et defined hy t = constant 
and v.- is the covariant derivative on the hypersurface Sj. From the determinant of the metric ([T|) one obtains 



3 



For general Hofava-Lifshitz-like gravity models, we do not require the full difFeomorphism-invariance, but only 
invariance under "foliation-preserving" diffeomorphisms: 

Sx'^C{t,x), St ^ fit). (4) 

Therefore, there are many invariants or covariant quantities made from the metric, in particular ^iji V^- Kjk^ ■ ■ ■ , 

^f^^ ■ ■ ■ ^^K,,, . . . , Ri^\ i?g), R<^1 Vf )i?f,L, • • • ; Vi^) . . . v(>gL, • • • , V, (n^'V.n^ - n-^V^n^, • • • , 
etc. Then the general consistent action composed of invariants that are constructed from such covariant quantities, 



ij ' ijkl^ i jklra^ ' *i «2 jklm^ 

could be a rather general action for the generalized Hofava-Lifshitz gravity. Note that one can also include the 
(cosmological) constant in the above action. Here it has been assumed that the action does not contain derivatives 
higher than the second order with respect to the time coordinate t. In the usual F{R) gravity, there appears the 
extra scalar mode, since the equations given by the variation over the metric tensor contain the fourth derivative. By 
assuming that the action does not contain derivatives higher than the second order with respect to the time coordinate 
t, we can avoid more extra modes in addition to the only one scalar mode which appears in the usual F(R) gravity. 
For example, if we consider the action containing the terms like 

(V'^V^)"+^ , (V''Vp)"V^KV,n"^-n''V,n^) , (6) 

the equations given by the variation over the metric tensor contain the fifth or higher derivatives (for a review of 
Hamiltonian structure of higher derivative modified gravity, see Q ) . If we define new fields recursively 

the equations can be rewritten so that only second derivatives appear. The scalar fields in ([7|), however, often become 
ghost fields that generate states of negative norm. Thus, we only consider actions of the form given by ([5]) in this 
paper. 

In the Hofava-Lifshitz-type gravity, we assume that N can only depend on the time coordinate t, which is called 
the projectability condition. The reason is that the Hofava-Lifshitz gravity does not have the full diffeomorphism- 
invariance, but is invariant only under the foliation-preserving diffeomorphisms UN depended on the spatial 
coordinates, we could not fix N to be unity {N = 1) by using the foliation-preserving diffeomorphisms. Moreover, there 
are strong reasons to suspect that the non-projectable version of the Hofava-Lifshitz gravity is generally inconsistent 
01 . Therefore we prefer to assume that N is projectable. 

In the FRW space-time with the flat spatial part and the non-trivial lapse N{t), 

3 

ds^ = -N{t)Mt^ + a{tf {'^^'f ' (8) 



i=l 



we find 



AT n^H 
Ko = irr^ ^.^^h. ^)o = HS\ other r(;=0. 



^.. = ^%, vf^^O, i?g,=0, V,{n^V.n^-n'^V.n^)^^^(^] , (9) 



where H = - \s the Hubble parameter. Then one gets 



R''' - - B^l - ^f^R^rn = = vifVg) . . . . . . . = , (10) 
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and since F must be a scalar under the spatial rotation, the action ([5]) reduces to 

Therefore, if we consider the FRW cosmology, the function F should depend on only two variables, ^ and 

3 d { a^H \ 
IiFN dt \ N )■ 

As a specific example of the above general theory, we may consider the following modified Hofava-Lifshitz F{R) 
gravity, whose action is given by 

= ^ / A''x^/^)NF{R) , R = K'^K.j - XK^ + 2^iW^ {n^^V .rf - n^V ,n^) - cf (gg') . (12) 

(3) . . ... f3) 

Here A and /i are constants and £^ is a function of the three-dimensional metric g^j and the covariant derivatives 

Vj-"^-* defined by this metric. Note that this action (IT2]) is more general than the one introduced in ref. ^] due to the 
presence of the last term in R. We normalize F{R) or redefine so that 

F'(0) = 1. (13) 

In 3, 

is chosen to be 

cf{g',f)^E^'g.,ME'^', (14) 
where Gijki is the "generalized De Witt metric" or "super-metric" ("metric of the space of metric"), 

gijkl ^ 1 ^gi3)tk g{3)jl ^ gi3)ilg(3}jk'^ _ Ag^^^^^^f^^'"' , (15) 

defined on the three-dimensional hypersurface T,t. E^^ can be defined by the so called detailed balance condition by 
using an action on the hypersurface Et 

^E^.^'-^, (16) 

and the inverse of Q^^'^'- is written as 

yijki - 2 {9ik 9ji + 9a 9jk J - ^9ij 9ki ' ^ ~ ' ^ ' 

The action W^[gj.; ] is assumed to be defined by the metric and the covariant derivatives on the hypersurface Ej. There 
is an anisotropy between space and time in the Hofava-Lifshitz gravity. In the ultraviolet (high energy) region, the 
time coordinate and the spatial coordinates are assumed to behave as 

x-^hx, t-^bH, z = 2,3, (18) 
under the scale transformation. In is explicitly given for the case z = 2, 

fd^x ^ _ 2A^) , (19) 

and for the case z = 3, 

f W3(r) , (20) 

where 

c.3(r) = Tr f r A dr + A r A r) ^ e^'^ (r^ld.ri^ + ^rr^r'^r^,) d'x . (21) 
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Here kw in ^T9\i is a coupling constant of dimension —1/2 and in ([20]) is a dimensionless coupling constant. A 
general E'^^ consist of all contributions to W up to the chosen value z. The original motivation for the detailed balance 
condition is its ability to simplify the quantum behaviour and renormalization properties of theories that respect it. 
Otherwise there is no a priori physical reason to restrict C^^^ to be defined by p^ . In the following we abandon the 

(3) 

detailed balance condition and consider Cj^ to have a more general form, since it is not always relevant even for the 
renomalizability problem. 

We now investigate the renormalizability and the unitarity of the model (|12p. For this purpose, by introducing an 
auxiliary field A, we rewrite the action (jl2p in the following form: 



= 2^ / d^xv/^iV {F'iA)iR -A)+ . (22) 

For simplicity, the following gauge condition is used: 

N = l, N' = 0. (23) 

Then one finds 

ro --io(3) p - io(3)''=o(3) r\ - r^^)' = io(3)./ f (3) , o) „„(3)\ 

other components of r(^p = , (24) 

and therefore 

(nn = (1,0, 0, 0) , K,, = ^glf , (n'^V.n'^ - n-^V.nn = (g^'^'' g^f) + \ (g^'= i^)" ■ (25) 



We define a new field by 



<^^ihiF'(A), (26) 



which can be algebraically solved as A = A((/?), so that 



The spatial metric is redefined as 



V = UnF'{A{^)) ^ F'iA{ip))=e'^. (27) 



9^^=e-''gf^. (28) 



Then the action has the following form: 



i + ^ - I) 9^^^^9l?^ + (? - T + f) + (^g ' - ^(-)} ■ 



Here 

4) (sg), ^) ^ (e-^g?) , l^(^) ^ A (^) F' (A (^))) - F (A (^)) . (30) 
If we insert ip = 1 into the action (I^Hl) . the standard Hofava-Lifshitz gravity emerges. On the other hand, if we choose 



M = A-i, (31) 



ip decouples with g^^^ . When the decoupling ([3T|) is assumed and 



A > ^ , (32) 



6 



Lp becomes canonical and the theory becomes unitary. In the case 

A - ^ , (33) 

the <y9^ term vanishes and therefore ip becomes non-dynamical, i.e. an auxiliary field. Eq. pi[) tells that /i = when 
A = 1/3. 

In order to clarify the renormalizabilty issue, we need to explicitly construct (^^i^) OS- ^ model 
corresponding to 2 = 2 in ()18p . which is still not renormalizable, we may propose 



(9f)=c.[^-I^+a(l^^)^ , (34) 



where C2 and a are constants. Since the action (p4|) induces the higher derivative terms to contribute to the propagators 
and therefore the propagators behave as 1/ |fc|* in the high energy region, the ultraviolet behavior is improved, although 
the theory still is not renormalizable. 

By the scale transformation ([^5]) . the curvatures are transformed as 



d(3) _ p(3) 



+ 2A(3)^ _ ig(3)Hv(3)^^(3)^) (35) 



Here R^f, R^^\ vf\ and A are the Ricci curvature, the scalar curvature, the covariant derivative, and the 
Laplacian defined by the metric g^j , respectively. Then if we consider the perturbation from the flat background. 



where Lp ^ due to ([13 



(3) 



(36) 



we find 



^ [d.d^hf - A/^i^)^}' + + 4a) {a.9^7.g) - A^^)''^} + Q + 4a) (A^.)^ 



Therefore if one chooses 



(37) 



(38) 



ip decouples with /ig\ Eq. ([37)) shows that the propagators of ip and /ig^ behave as 1/ in the high energy region, 
so that the ultraviolet behavior is improved. 

Similarly, a model corresponding to z = 3 in ([TSl . which could be power-counting renormalizable, can be obtained 
by choosing 



(^g)) = C3 (^V(3)'=i?(3)^^ vg)i?g) + iv(3)'^-i?(3)vg)i?(3)) . (39) 



For the z = 3 model, the dimension of ip vanishes and therefore all the interactions in (di^^ ^ V') ^^^d V{ip) in ([29[) 

become power-counting renormalizable. The propagators of ip and ft.g'' behave as 1/ |A;|^ in the high energy region, so 
that the ultraviolet behavior is improved to be renormalizable. 

We have shown that by requiring ([?T|) and (1551) . the scalar field decouples with the gravity modes in the Einstein 
frame. The decoupling itself is not directly related with the renormalizability, but the decoupling makes it much 
easier to discuss the renormalizability of the model. The choice p9p for gives the renormalizable model. The 
renormalizability does not essentially depend on the functional form of F{R). 



7 



III. HAMILTONIAN ANALYSIS OF THE GENERAL ACTION IN THE FRW SPACE- TIME 



Let us analyze the proposed general action ([Tl]) of the FRW space-time ([8]) with the flat spatial part and the 
non-trivial lapse N — N{t). Introducing four auxiliary variables a,A,l3,B enables us to write the action PT|) as 



5gHL 



a A 



H 

N 



/3[B 



3 d /a^H 



a^N dt\N 



FiA,B) 



The variations of the action (|40p with respect to a and /? yield 

3 d 



A 



H 

N 



and B 



a^N dt \ N 



(40) 



(41) 



respectively. Integration by parts permits the removal of the second-order time derivative of a and the time derivative 
of iV, assuming the boundary terms vanish, but with the price that /3 becomes a dynamical variable. Thus the action 
(Uni) can be written as 



5, 



gHL 



d^xVgi^N 



(42) 



The action is equivalent to (HO]) and consequently to the original action ((TT|) . The advantage of the action (H^ 
over is the simpler dependence on the variables a and N , which will be crucially important in the following 
Hamiltonian analysis. 

For the Hamiltonian analysis of constrained systems and their quantization we refer to the monographs lol- 13|. 

In the Hamiltonian formalism the generalized coordinates gfj ■, N, a, A, (3 and B of the action (|42l) have the 
canonically conjugated momenta tt'-', tt^v, tt^, tt^, tt^ and ttb, respectively. We consider N to be projectable, N = N(t), 
and therefore also the momentum tt^v — T^N(t) is constant on the hypersurface St for each t. The Poisson brackets 
are postulated in the form (equal time t is understood) 



{a{x),TTa{y)} = 6{x-y), {A{x),'KAiy)} = S{x - y) , 
{(3{x),TTfs{y)} = S{x - y) , {B{x),TTBiy)} = S{x - y) , 



(43) 



with all the other Poisson brackets vanishing. We are considering the FRW metric ^ with the flat spatial part 
gf^ = d^bij and therefore the Poisson bracket for the scale factor a and the momenta conjugate to the 3-metric takes 
the form 



d3y{a,^*^(y)} 



2^ 



(44) 



Let us find the momenta and the primary constraints. The action (j42p does not depend on the time derivative of 
iV, a, A or B. Thus we have the primary constraints 



<I>i=7rjv«0, $2(a^) = 7r„(a;) « 0, $3(0;) = 7rA(a;) w , <^ ^(x) = -k b(x) k. Q . 

(3) 



The momenta conjugated to /3 and g\i are 




(45) 

(46) 
(47) 



respectively. The "velocities" /3 and can be solved in terms of the canonical variables, so there are no more 
primary constraints. 



8 



Then we define the Hamiltonian 



H 



(^'^g^f + ^^/j) -L= f d'^xNU , (48) 



where the Lagrangian L is defined by p2)) . S'gHL = / dti, and the so-called Hamiltonian constraint is found to be 



^ = Y (^^ I] ^" + - (aA + pB + F{A, B)) . (49) 



The primary constraints (|45p can be included into the Hamiltonian (|48p by using the Lagrange multipliers Xk, k — 
1, 2, 3, 4. We define the total Hamiltonian by 

4 p 

Ht = H + Ai$i + ^ / d'xXn{x)<^n{x) . (50) 

ri=2 

Note that there is no space integral over the product Ai<i>i = Aittjv, since these variables depend only on the time 
coordinate t. 

The consistency of the system requires that every constraint has to be preserved under time evolution. Since the 
Poisson brackets of the primary constraints ps)) are zero, 

{$,.,$,1 = 0, fc,Ze {1,2,3,4}, (51) 

the time evolution of the primary constraints is determined by the Hamiltonian H alone 

$fc ^ Ht} = i/} , fc = 1, 2, 3, 4 . (52) 

Thus the following time derivatives of the primary constraints have to vanish: 



$1 


= T^N 


= {7rN,H} 


$2 


= T^a 




$3 


= T^A 




$4 







,3 



d^xH 

dF{A,B] 



3 



a 



OA 
dF{A,B) 



dB 



(53) 



Since none of these expressions ([53)) vanish due to the primary constraints p5|) . we must impose the secondary 
constraints: 

$0 = / d?xH « , 



Here the position argument x has been omitted in the right-hand side of the local constraints. Note that neither 
N or a can be constrained to vanish, since they are the essential physical quantities in this theory. Here the actual 
Hamiltonian constraint $o is global due to the projectability condition, N = N{t). Note that the Hamiltonian (|48|) 
is simply this constraint multiplied by N, i.e. 

H = N<^o . (55) 
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Also the secondary constraints ([54)) have to be preserved under time evolution. The time evolution of the secondary 
constraints is 

- iJr} = iV{$,„, $o} + ^ / d^y A„(y){$,„, $„(y)} , m = 0, 5, 6, 7 , (56) 

n=2 •' 

where we have used (|55p and the fact that none of the constraints = 1, 2, . . . , 7, depend on the lapse TV, and 
that the secondary constraints (|54p do not depend on tt^v- For the global Hamiltonian constraint $o we find the 
following Poisson brackets with the primary constraints (l45t 



{$o,*2(a;)}--3'5(a;), {$o, ^sN} = -a'$6(a;) , ^ ^a^^j{x) , (57) 



which all vanish due to the other secondary constraints. Thus, according to ([56]) and ([57]), the Hamiltonian constraint 
$0 is preserved under time evolution, $0 ~ 0. For the secondary constraint $5 we obtain the non-vanishing Poisson 
brackets with the primary constraints (1451) and the Hamiltonian constraint $0 

{^^{x),^:i{y)} = aH{x-y), {$5, $o(aj)} = + 3^^A . (58) 
For the next secondary constraint $g we obtain the non-vanishing Poisson brackets 

{$6(a;),*2(y)} = 5{x-v), 

{Mx),^,iy)} = ^^£^^5{x-y). (59) 
For the last secondary constraint $7 we obtain the non- vanishing Poisson brackets 

{<f,ix),^,iy)} = ^^l^^6{x-y), 

{$7(^),<i>4(y)} = ^^^^5{x-y), 

{$7(0;), $0} - \ (^Jy''+<^ ■ 

Inserting all these Poisson brackets into ([55[) gives the tertiary constraints 



(60) 



$5 = N (^-^ + STrpA^ +X3a^ ^0. (61) 

$6 = X2 + X3 +A4 g^g^ «0. (62) 

• TV /2 A \ ^ d^F(A, B) , d'^F{A, B) „ 

We assume that all the second partial derivatives of F{A, B) do not vanish.^ In this case the equations ([6n) - ([63[) are 
restrictions on the Lagrange multipliers, constituting an inhomogeneous linear equation for the unknown multipliers 



^ This is the case for example in the modified Hoi'ava-Lifshitz gravity model F{R) on R + bR^ discussed in that corresponds to 
F{A, B) = F((3 - 9A)A2 -I- 2^^B) oc 6(3 - 9\f A'^ + 2h^l{Z - <i\)A^B + ibfi^B'^ + (3 - 9X)A^ + 2^lB , 
so that we would have 

d^F(A,B) ,.00, , . d'^F(A,B) , d^F(A,B) , 

oc 126(3 - 9\fA^ + ibfj.(3 ~ 9X)B + 2(3 - 9A) , oc Ab^J.{3 ~ <i\)A , oc Sbfj.^ . 

oA^ oAoB aB^ 
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Xi,i = 2,3,4. Since the homogeneous part of this equation has only the null solution A2 = A3 = A4 = 0, the most 
general solution is the solution of the inhomogeneous equation: 



A. 



N 



A, = 



N 

3 

i=l 

N 



B 



_ %TTiiA \ d^F{A, B) 
J dAdB 



Nm=-\B 



a+\B 

d-KpA 



d^F{A,B) fd^F{A,B) 



dAdB 



N 

Xa = Nua = — 



9 3 



9^pA\ d^F{A,B) 
"+'^"^J dAdB 



dB"^ 



d^F{A,B) 
dB^ 



(64) 



The multiplier Ai is arbitrary, as is the non-dynamical variable N that also is a multiplier in the Hamiltonian ([ST 
with ([55)1 . 

The total Hamiltonian (|50p can be written as a sum of two first-class constraints multiplied by the two arbitrary 
time-dependent multipliers N and Ai: 



Ht = NHq + Ai$i 
where we have defined the first-class Hamiltonian constraint by 



4 r 

+ H / d^a;M„(a;)$„(a;), 

n=2 •' 



(65) 



(66) 



with the fields u„ [n — 2,3,4) given by It is easy to see that $1 = ttjv is first-class, since it clearly has a 

vanishing Poisson bracket with every constraint. From (1521) and (j56p we see that the sum of constraints is first- 
class by construction. Note that (j66p is a combination of secondary and primary constraints. Usually a secondary 
first-class constraint would require us to define an extended Hamiltonian where the constraint would be added with an 
additional arbitrary multiplier. In this case, however, that would only lead to a redefinition of the multiplier N , and 
such a change W — )■ TV -|- an arbitrary function of time' does not bring anything new to the description. As always, 
the first-class constraints are associated with the gauge symmetries of the system [l6| . The first-class constraints i/o 
and $1 generate the (gauge) transformations that do not change the physical state of the system. 

The constraints Xk = ($2, 'I'a, ^'5, ^"67 'I't) form the set of second-class constraints of the system.^ For details 
on the classification and representation of second-class constraints, one can see 
second-class constraints define the matrix: 



17j |. The Poisson brackets of the 



Cki{x,y] 



where 



Cki (x) 



= {Xk{x 


),xi{y)} = 


Cki{x)S{x 




/ 





-1 


^ \ 





-a^ 


-Fa2 - 


Fab 








-Fab - 


-Fb2 











1 
3 


1 Fa2 


Fab 








\0 Fab 


Fb2 -i 





J 



(67) 



(68) 



and we denote 



d'^FjA, B) 
dA^ 



Fab 



^ d^F{A,B) 
^ dAdB 



_ d '^FjA, B) 
dB^ 



(69) 



^ The index of Xk runs over fc = 1, 2, . . . , 6, so that Xk = 'I'fc+l- 
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This matrix has the inverse 



C^\x,y) = C^\x)5{x-y), (70) 





Fae 



3^F 2 

F^2F^2-FX, 
a^F„2 



Far 

F„2 



Fab 
3a^F„ 



1 

3a-iF„ 



3a-' Fg 
1 

3a^Fg2 









1 

Fab 







- F^2F^2 -, _Fab \ 

1 _ _ 













(71) 



which satisfies 



J d^^ Cki{x, z)&^{z, y) = Cki{x)&'^{x)5{x - y) = ^™^(a: - y) . 



(72) 



Now it is possible to impose the second-class constraints Xk by replacing the Poisson bracket with the Dirac bracket. 
For any two functions or functionals / and h of the canonical variables, the Dirac bracket is defined by 

{/(a;),My)}DB ={/(a;),My)}-/dW^'{/(a;),Xfc(^)}C^-'(^,^'){x;(^'),My)}- (73) 

The Dirac bracket takes fully into account how the second-class constraints impose relations between the canonical 
variables. Therefore it enables us to set these constraints to vanish strongly, Xk{x) — 0. So we have the identities 



TJ'a = = TTS = , A = 



3a3 



and 



dF{A,B) 



OA 



, /3 



dB 



(74) 



(75) 



When the function F is known, from ([75|) we can solve the variable B in terms of /? and ^ = 



Then a can be solved: 



dA 



(76) 



(77) 



Introducing these strong constraints into the Hamiltonian gives 

3 



3a 



E 

i=l 



TT —a 



(78) 



The first-class Hamiltonian ([M)) reduces to Hq = $o and the total Hamiltonian becomes 

Ht = N<i>o + Xi^i = N J d^x'H + XiTTN ■ (79) 

The canonical variables are iV, ttjv, (7,||-',7r*-', /3,7r^. In other words a,A,B and their conjugated momenta have been 
eliminated. 

In order to obtain the equations of motion, 



/ = {/, ^^tIdb = N{f, $o}db + Ai{/, ^atIdb : 



(80) 
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for the canonical variables we have to work out all the Dirac brackets (1731) between the variables. We find that 



the Dirac bracket (|73|) reduces to the Poisson bracket (|43l) for all the canonical variables iV, ttjv, '/ll'', tt*^ , /3, tt^, and 
consequently for any functions of these variables. In the first pair of variables, A'^ is quite arbitrary and ttn does not 
evolve due to the equations of motion: 



N = {N, Ht}bb = Ai{iV, ttn} = Ai , 

TTjv = {nN,HT}DB = {ttat, N} J d^xn = - y d^xn 



0, 



where as before Ai is an abitrary function of time. For the spatial metric we get 

9i 



•(3) r (3) Tj 1 2Ntti3 



3a 



where g^-^-* — 2aaSij. Solving for a gives 



/ dtNiTf). 

J to 



(81) 



(82) 



(83) 



Hence we need in order to get a{t). This reveals that tt^ does not depend on the spatial coordinate x, because 
both a and N depend only on the time coordinate t. For the conjugated momenta we obtain the equations of motion 



3a3 



k=l 



kk 



^ Ub + f{a,b) 



TTff dF{A,B) 
2a2 dA 



(84) 



where the arguments of B are omitted for brevity, B = B{/3,A) = B (^/3, ^) , as will be in the next equation. For the 
variable /3 we obtain the equation of motion 



dF{A,B) 



dA 



3o^ 



For its conjugated momentum tt^ we obtain the equation of motion 

7T/5 = {7r^,i7T}DB = Na'B ^) 



(85) 



(86) 



Further progress in the study of dynamics practically requires one to specify the form of the function F ^ and then 
solve ([75]) from ([75]). We can conclude that when the second partial derivatives of the function F are non-zero, 
the proposed general action (1111) defines a consistent constrained theory. 

Let us then briefiy consider the cases when some of the second partial derivatives (j69p of the function F are zero. 
In such cases the tertiary constraints (pT|) -(|M | are no longer mere restrictions on the Lagrange multipliers, but in 
addition impose constraints on the canonical variables. As an example we consider the case when F^i — and 
-F42 7^ 0,-Fab 7^ 0. Then we obtain the tertiary constraint 



•iTxpA\ 



Fab^O; 



(87) 



and solve two of the Lagrange multipliers, say A3 and A4: 



A, = AT 



B 37r^^^ 
3" ^ ~ 



A4 = -- 



Fab 



X2+N 



B 3TrpA\ 



Fa2 



where the third multiplier A2 is arbitrary. The consistency condition w of the tertiary constraint imposes 
a quartic constraint on the canonical variables, because $§ turns out to be independent of the Lagrange multiplier 
A2 and non-vanishing due to the constraints established so far. Further constraints may follow from the consistency 
condition of the quartic constraint. This has to be checked explicitly after choosing the form of the function F. These 
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additional constraints are a serious threat to the viabihty and consistency of the action, since they may delete the 
physical degrees of freedom. In case we also have F^b = 0, we would solve the Lagrange multipliers A2,A3 from 
(pT|) - ([S5)) and obtain a tertiary constraint that restricts the field /3 to be a constant, /3 « 0. Thus in the latter case 
we should not have introduced the auxiliary fields /3 and B in the first place, since F in the action is already linear 
in its second argument. We do not discuss the case ,^42 = 0, because it appears to have very little if any practical 
application. 

As a specific example of the above general theory, one can consider the FRW cosmology in the modified Hofava- 
Lifshitz F{R) gravity studied in and its further generalization considered in the present work, as action (|12p. 
However, this analysis can be used to study FRW cosmology in any theory with an action of the general form ([S|). 
Moreover, the methods presented in this section can be used to analyze any action of the form ([5]) in a general way, 
without assuming any particular space-time. The proposed modified Hofava-Lifshitz F{R) gravity will be studied in 
the next section. 



IV. HAMILTONIAN ANALYSIS OF THE F{R) GRAVITY 



Let us then consider the Hamiltonian analysis of the proposed action for the modified Hofava-Lifshitz F{R) 
gravity. The analysis is similar with the Hamiltonian analysis presented in ref . [5|] , where a special case of this theory 
given by the choice was proposed (see also the analysis of ref. IJ]). This special case with the further restriction 
to the parameter value fi — has been proposed and analyzed in ref. [l5| . In this section we generalize the analysis 
of ref. 3- 

By introducing two auxiliary fields A and B we can write the action (|12p into a form that is linear in R: 



S 



FiR.) 



5(3)iv B{R-A)+F{A) 



(89) 



Then we can write R as 



Introducing (|90p into (j89p and performing integrations by parts yields the action 

Sp^j,) = I dtd^x^{N [b (K^.g^^^'Ku - Cf (sg)) - a) + F(A) 



13- N'-d^B) - 2iiNg 



,(3)»jy(3)y(3)^ 



}■ 



(90) 



(91) 



where the integral is taken over the union hi of the t — constant hypersurfaces Sf with t over some interval in M. We 
assume that the boundary integrals on dlA and dT,t vanish. The difference compared to the action studied in ref. Q 
is that the potential part (^g"* ) may have any form that satisfies the correct scaling property under (fTS]) . In other 
words it is not necessarily defined by and the detailed balance condition p^ . However, due to the projectability 



condition N ~ ^(1) the specific form of the {gl^ ) has very little effect on our analysis. Indeed the analysis of 
ref. [sl is translated to the present more general case by making the replacement p^ from rhs to Ihs. Therefore we 
only present the main points of the generalized analysis. 

In the Hamiltonian formalism the field variables gij^ N, Ni, A and B have the canonically conjugated momenta 
TT*-', TTjv, tt', tta and tt^, respectively. For the spatial metric and the field B we have the momenta 



SS 



F(R) 



Sgij 
SS 



BQ^^'^'Ku - ^5^'^'-' ( B - N'd,B 



F(R) 



SB 



(92) 
(93) 



We assume /i ^ so that the momentum ((93|) does not vanish. Because the action does not depend on the time 
derivative of N, or A, the rest of the momenta form the set of primary constraints: 



TTJV 



0, tt'{x)kO, TrA{x)^0. 



(94) 
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Due to the projectability condition also the momentum is constant on Et for each t, ttn = 7rjv(i)- Then the Hamiltonian 
is calculated 



H 



(95) 



where the so-called Hamiltonian constraint and the momentum constraints are 



;(3) 

7(3), 
j 



Hn — 



B (Cf ) +a)- F{A) + 2m<?(3)'^ Vf Vf 

respectively. We define the total Hamiltonian by 

Ht = H + Xnttn + j d^x (AiTT^ + Xa^a) , 



B 



1 -3A 



-BTxi 



(96) 



(97) 



where the primary constraints (|94p are multiplied by the Lagrange multipliers Aat, Ai, \a- 

The primary constraints (|94p have to be preserved under time evolution of the system. Therefore we impose the 
secondary constraints: 



<i>o = J d^xHo^O, $l;(a;)EE-H'(a;) wO, ^ a{x) = B (x) ~ F' {A{x)) . 



(98) 



Here the Hamiltonian constraint $0 is global and the other two, the momentum constraint $5(2;) and the constraint 
$^(a;), are local. It is convenient to introduce a globalised version of the momentum constraints: 



(99) 



where i = 1, 2, 3 are arbitrary smearing functions. It can be shown that the momentum constraints $5(^1) generate 

(3) 

the spatial diffeomorphisms for the canonical variables -B, ttb , glj , tt*^ , and consequently for any function or functional 
constructed from these variables, and treats the variables A, tta as constants. 

The consistency of the system requires that also the secondary constraints $0^ ^s{£,i) and ^a{x) have to be preserved 
under time evolution defined by the total Hamiltonian (197^ . which can be written in terms of the constraints as 



Ht = N^o + ^s{N^) + Xn^n + J d^x (AiTr* + Xa^a) ■ 
The Poisson brackets for the constraints $0 and ^s{£,i) are 

For the constraints tta and ^a{x) the Poisson brackets that do not vanish strongly are: 

{7:a{x), $0} = - V7^*A(a;) « , {7rA{x),<^>Aiy)} - F"{A{x))d{x - y) 



(100) 



(101) 



{<^o,'^a{x)} = 



1 -3A 
2/i 



BtTb] , {<^>s{^^),<^A{x)} = -e^^B. 



(102) 



Since F"{A) — would essentially reproduce the original projectable Hofava-Lifshitz gravity, we assume that ^"(^4) 7^ 
0. The constraint ^a{x) can be made consistent by fixing the Lagrange multiplier A^: 



Xa 



1 



F"{A) 



N'diB 



N 



(103) 



15 



Now all the constraints of the system are consistent under dynamics. 

According to the Poisson brackets (|10ip - (ll02[) between the constraints, we can set the second-class constraints 7r^(a;) 
and <l>^(a;) to vanish strongly, and as a result turn the Hamiltonian constraint $o and the momentum constraint ^siS,i) 
into first-class constraints, by replacing the Poisson bracket with the Dirac bracket. It turns out that the the Dirac 
bracket reduces to the Poisson bracket for any functions of the canonical variables. Assuming we can solve the 
constraint $^(a;) = 0, i.e. B = F'{A), for A = A{B), where A is the inverse of the function F' , we can eliminate the 

(3) 

variables A and tta- Thus the final variables of the system are glj , tt'-', B, ttb- The lapse N and the shift vector Ni, 
together with and A^, are non-dynamical multipliers. Finally the total Hamiltonian is the sum of the first-class 
constraints 

Ht = iV$o + «'s(A^i) + Aatttat + J d^xX,T:' . (104) 

We conclude that also the proposed action of the more general modified Hofava-Lifshitz F{R) gravity defines 
a consistent constrained theory when the projectability condition is postulated. For additional details and discussion 
on the analysis see ref. j5|]. 



V. HAMILTONIAN ANALYSIS OF THE F{R) GRAVITY IN FIXED GAUGE 



Let us then analyze the action (IT^ when the gauge is fixed by and we obtain the action (P^ . First we find 
the momenta canonically conjugated to the generalized coordinates g^-J'' and (p of the action (I^Hl) . For the fields (p and 
9i? ^® ^^'^ momenta 



5S 



F{R} 



and 



Sip 4/t2 



(-(1 - 3A + 3M)ff(')'^^g,|f + 3(1 - 3A + 6m)(p) 



respectively. 

In the following analysis we will first assume 



l-3A + 3Ai7^0, l-3A + 6Ai7^0, fi ^ , 



(105) 



(106) 



(107) 



so that the kinetic term for tp docs not vanish. Later we will consider the special cases where these conditions do not 
hold. First we solve (p from (|105p . 



4k-' 



=7r^ -f (1 - 3A + 3/i)5(')'^'5lf 
3(1 - 3\ + 6fi) 



and introduce it into (|106l) 



5(3) 



4k2 



9ki 



1 



3 l-3A + 6Ai 



Then we find the velocities in terms of the coordinates and momenta. First we contract p09p by g\, and solve for 



l-3A + 3/i (3^,, 
3(l-3A-f6M)^ ^' 

(3) 



(108) 



(109) 



4k^ 



'g 



(3) 



1 - 3A -t- 6/i_(3) 



1 - 3A -t- S^i 



(110) 



"^'^ 9a*2 

This is inserted back into (|109p as well as into (|108p , which enables us to obtain the velocities in terms of the canonical 
variables: 

4^2 



ip : 



'9' 



3A-1 1 -3A-h3M_(3). 



(3) V 27/i2 ^ 



^(3) 

91, 



7(3) 



-A3) 



1-3A + 6M lV(3)-fei l-3A + 3Ai 



3 ] 9ki ^ 



(111) 
(112) 
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Thus there are no primary constraints. It is expected that there are no first-class constraints, since the gauge has been 
completely fixed by setting A*" = 1,7V* = 0. Due to the non- vanishing kinetic terms of (/?, there are no second-class 
constraints either. The Hamiltonian is defined by 



H = (Fx 



(if'' 



^(3) 



(113) 



After a lengthy algebra exercise we find 



H 





4^2 











1 



1 - 3A + Ba* _(3) 



3A- 1 



54/i' 



2 V 



'9 



(3) 



2k2 



(3) 



(114) 



In fact wc find that the Hamiltonian (|114p is correct for any parameters A and fi as long as it is defined, i.e. when 
/i 7^ 0. This can be seen by considering the two cases when the kinetic cross-term vanishes {(p and g^^ decouple), 
1 — 3X + ~ 0, fi 0, and when the kinetic (p'^ term vanishes, 1 — 3A + 6/i = 0, /i 7^ 0, separately. Even the formulas 
(|llip - pi2p hold in these cases, thought the details of their calculation are quite diferent. Note that the vanishing of 
the both kinetic terms of ip implies fi — and A = 1/3. 

The Poisson bracket is postulated by (equal time t is understood) 



{-g^f{x),7r'^\y)} = - {S'^S^ + 6iS^) S{x - y) , {p{x),7r^iy)} = S{x y) 



(115) 



with all the other Poisson brackets vanishing. Now we can work out the Hamiltonian equations of motion. Because 
there are no constraints, the Hamiltonian (|114l) defines the dynamics of any function or functional / of the variables 



,(3) 



, tt'^ , (p, TT^ by: 



f^{f,H} 



(116) 



For the generalized coordinates glj and we obtain the equations of motion (|112p and (jllll) respectively. For the 
momenta tt*-' we obtain 



4^2 



1 



.(3) (3) -feZ-mn 



a f_(3)-kl\' 

4 i^^' " ) 



-9ki 



\^) -ik -jl I ;=:i7-(3)-fcZ I h-ij ^ 

1/tt tt-' + air gl.iTT ^-Ott-'tt, 



-p) / 1 



2k2 



±o(3)*J 

2^ 



L 



(3) 



(g(f,^)-y(v.) 



dg, 



(3) 



(117) 



where we have introduced the constants: 

1 - 3A -f e^i 1 



27/i2 



b = 



1 - 3A + 3/i 
2V 



The equation of motion for tt^, is 



3A- 1 
27/i2 



2^2 



dp 



3A{ip)e^^ , 



(118) 



(119) 



where the last term is obtained from the derivative of V{(p) from ([50]) 

dy(^) ^ ^(^)dAM^„(^(^)) ^ 3Aip,)F'iA{^)) ^ 3A{^)e'- 



dip 



dip 



(120) 
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Here we have also used the definition ([26)) of (p and ([27)1 of A{lp) to calculate 

In particular, the equations of motion (I117P for the momenta tt*^ are pretty complex. However, as always, they are 
first order differential equations. 

The cases with fj, = are less interesting and we only consider them briefly. When A = 1/3 the field is non- 
dynamical and hence one has the primary constraint tt^ « 0. The momenta conjugate to g^^- are given by 

= ^ {9'''"'-9'''''-9^^ - 1-9^'^'' 9''^'' -9^^^) ■ (122) 

It has zero trace g|f^7f*^ = and can be trivially solved for g^f = 1" glfc''g!F^^'- When A ^ 1/3, one is forced to 
impose the constraint tt^ — —g^^'^Tt''-' ~ that again leads to (jl22p and makes ip non-dynamical. 



VI. FRW COSMOLOGY IN POWER-LIKE MODELS: COSMIC ACCELERATION AND FUTURE 

SINGULARITIES 

We now consider the FRW cosmology of the action (fT2|) . In the spatially-flat FRW space-time ([8|), since the spatial 
curvature vanishes, = R^^^ — 0, there is no contribution from >C^\ as it vanishes according to or ([M)) . In 

other words, the choice of Cj^ in (l34t or p9l) gives the same FRW cosmology. Of course, this situation changes when 
one considers black holes or other solutions with non-trivial dependence on the spatial coordinates. 

Let us first review the spatially- flat FRW equations obtained in ref. [5]. Varying the action (fT2)) with respect to 
and setting iV = 1 one obtains: 



= F - 2 (1 - 3A + 3/^) (ij + 3H^^ F' {r^ -2{l-3X)H 



dF 



d^F' R 



dt 



2^- 



dt2 



(123) 



where F' denotes the derivative of F with respect to its argument. Here, the matter contribution (the pressure p) is 
included. On the other hand, the variation over N gives the global constraint: 



= / d' 



dF' [R 



dt 



(124) 



after setting = 1. Here p is the energy density of matter and we have set again = 1. It is important to stress 
that, because of the projectability condition N — N{t), the above equation is a global constraint. If the standard 
conservation law is used, 



0^p + 3Hip + p) , 



(125) 



Eq. (|123l) can be integrated to give:' 



= F i? -6U1-3A 



|(l-3A + 3At)i?^+M-ff}i^' (-R) +6pH 



dF' 



'(«) 



dt 



C 



(126) 



^ Note that, as already shown in [Tsl l for the standard case, the parameter A has a crucial role in the relation of Hoi'ava-Lifshitz type 
theories. In fact, from the second equation above one realizes that this solution is physical only if 1 — 3A -|- 3/i > 0. It should also be 
pointed out that in Hofava-Lifshitz gravity the role of standard matter and its conservation properties are not well understood yet. We 
will proceed with our discussion supposing that it is possible to couple matter and gravity in the same way in which one docs in GR. 
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Here C is the integration constant and can be set to zero. In [19[, however, it has been claimed that C does not 
necessarily need to vanish in a local region, since (|124p needs to be satisfied only in the whole universe. In this sense 
in a limited region, one can have C > and the Ca~^ term in (|126p can be regarded as dark matter. 

Note that Eq. (|126p corresponds to the first FRW equation and p23p to the second one. Specifically, if we choose 
A = /i = 1 and C 0, Eq. (fHS)) reduces to 

/ -x /~\ dF'(R\ 

^ F (r) -6(h^ + H] F' (r) +6H ^ ^ - 



dt 

(127) 

f(r] -6 (h^ + h) F' (r) + 36 (aH^H + h] F" (r) -p, (128) 



which is identical to the corresponding equation in the standard F{R) gravity (see Eq. (2) in 20] where a reconstruction 
of the theory has been made) . In the following we will explore the properties of the equations (|123p and (|126p , especially 
looking for solutions that represent accelerated expansion. Solutions of this type are very important because they 
represent the key evolutionary phases of the universe, namely the inflationary era and the dark energy era. The 
connection with dark energy is particularly important to understand, if the Newtonian nature of the quantum theory 
of gravitation, implicit in Hofava-Lifschitz gravity, is the direct cause of cosmic acceleration and, as a consequence, of 
dark energy. 



A. de Sitter cosmology 

Let us investigate the properties of the de Sitter solutions in this class of theories. This issue was considered for 
the first time in [5| , but in the following a more general treatment is given. These solutions are of great importance 
in cosmology because they have the potential to describe both inflationary phase(s) as well as dark energy era(s). In 
standard F{R) gravity it has been proven that it is possible to construct a viable model unifying inflation and late 



time acceleration in the form of double or multiple de Sitter solution j2l|-|24 
In vacuum {pm,Pm — 0) and substituting the de Sitter metric 



3 

> 2 



ds^ = -dt^ + exp (7t) ^ (dx') , (129) 
1=1 

the equations (|123p and (I126P reduce to the single equation 

F + 672(3A-3Ai-l)F' = 0. (130) 

In Table HI we show the values of the time constant 7 of the do Sitter metric of some popular F{R) theories and their 
Hofava-Lifshitz versions. 

It is interesting to note that, in contrast to what happens in standard F{R) gravity, the quadratic function F = i?"* 
is not degenerate for m = 2, but only for 

^ (131) 



1 - 3A + S^i ' 



i.e. it depends on the Lorentz- violation parameters. 

It is also interesting to note that in general the solution of (I123P and (|126p is not unique. Thus a given theory can 
have multiple dc Sitter solutions. This means that also in this case the cosmologies of these theories can admit both 
inflation and dark energy phases. However, since the Hofava-Lifshitz parameters are in principle only present in the 
coefficients of the equation (|130p and the number of solutions of (|130p is determined by the powers of R appearing 
in F, two corresponding theories will have in general the same number of de Sitter solutions. Obvious exceptions are 
the case 3A — 3/i — 1 = for which the Eq. (|130l) becomes F — and the case 3A — 3/i — 1 — 5(7)- 

In the first case, we see that the structure of the cosmological equations is essentially changed. In particular, the 
equation (|123p is modified and the constraint (|126p looses the linear term: 



dF'{R] 

= F (R) -6nHF' [R] +6pH - P (132) 
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(we have considered C = 0). Consequently, the equation (|130p becomes 

= F , (133) 

which is never obtained in the standard F{R) gravity. 

In the second case, instead, the choice of the function g can radically change the number and type of solutions in 
these theories. In this sense, the solution space for Hofava-Lifshitz F{R) gravity can be considered bigger than the 
one of its standard counterpart. Such a fact will be even more apparent in the case of the FRW-type solutions that 
will be examined in the next section. 



TABLE I: Some of the values of the time constant of de Sitter backgrounds for standard F{R) gravity models and their Hofava- 
Lifshitz counterparts. When writing the form of the function F, the Ricci scalar of both types of theories is indicated by x. 
For the more complex forms of F{R) an implicit equation has to be solved for the time parameter 7 in order to find its values. 



Function F Standard Case Hofava-Lifshitz case 



^ -rx ^1 _ n « _ n 

^ + X+ + -C + 67"+X = -3D + 67"(3A-3m + 1)+X = 

A = 2^<'"+"'3™+"(-m + n + 2)7^'™+"' + 2^™3™(2 - m)C7^'" + 2^"3"(n + 2)x7^" + 2{x 
B = 22('"+n)3m+n^2(m+,i) ^^-^ _ _ i) _ 3„A + 3n/i + n + 2) + 

-|-2^'"3"57^'" (m (3A - 3^ - 1) -I- 2) -|- 2^"3"x7^" (-3nA + 3nn + n + 2) + 2Cx 

p _ x(-6(-2re-3)-,25_i)_8(2n + l)(c, + l)725x 

^ _ 6(2re+l)(Q + l)T^5x(3A-3n-l)+x(-67=^5(n(6A-6n-2) + 3(A-M-l))-l) 
(l2-,2<;-l)(„(l272;_i)2^ + l^„^1^2 



B. Power law solutions and reconstruction technique 



In addition to the de Sitter solution described above one can also look for accelerated expansion phases in the 
form of power law solutions. These solution can have a double value as Friedmannian cosmologies, if the exponent of 
the power law is in the interval ]0, 1[, and they can realize the so-called "power law inflation", or a "power law dark 
energy" , if the exponent is bigger than one. 

If we look for the presence of Friedmann solutions of (|123p and ()126p . we realize quickly that, as in F{R) gravity, 
there is little chance to find power law solutions, unless one considers a function F of trivial form. However, due to 
the additional parameters, the set of solutions of this type is bigger in the Hofava-Lifshitz case than in the standard 
one. For example, in the simple case F{R) = R + xR™ we find that, in the presence of a barotropic fiuid {p ~ wp), 
the spatially-flat solution 

i2/3(l + lo) ,-2 /T3A\ 

a = aot ' ^ ^ ' , p = Pot , (134) 

satisfies dm)) and ([726]) if 

{w^-l)j{S\~l) 4(3A-1) 

■■ - and po = — . (135) 



^ 2w(3(u;-|-l)7-u; + l) ' 3(^-^1)2^ 

This corresponds to the standard Friedmann solution. It is well known that in standard F(R) theories, the case 
F{R) = R + xR"^ possesses only power law solutions of the type a = a^t or a = aot^/"^ (see e.g. 25|). 
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In order to facilitate the analysis in the next sections, we consider also some solutions for this model in the case 
of very small and very large scalar curvature. In the first case, the theory reduces itself to GR plus a cosmological 
constant and its solutions are approximated by the Friedmann ones. In the case of high curvature, instead, the theory 
reduces to F{R) w i?". Such a theory possesses three exact solutions. The first two 

a = aof^ , p = pot^^ , 7 = . , (136) 

3(1 + wj 

_ X [3/i(m - l)(2m(w + 2) - w - 1) - m{2m - 1) (3A - 1)] Mm^ (-3A + 6^* + 1) - 12m(w + l)p 
^ «;2 (to [3A - 6^ - 1) + 3{w + l)p] V (w + 1)2 

and 

9 2(m- 1)(2to- 1)(/ 

" = ^= (2^-l)(3A l)-6(r-l)M ' ^° = °' ^'''^ 

correspond to the solutions in the standard F{R) case. A third solution, that is valid only for m > 1, is 

a = aof^ , p = pot~'^ , 7 = — — - , Po = , (138) 

— 3A + D/i + 1 

which is characteristic of Hofava-Lifshitz gravity and does not depend on m. In the analysis of the singularities of 
this simple model, we will refer to these solutions. 

Note that as it often happens [i^ in theories of this type, the value of po can be negative (or even undefined) for 
certain combinations of variables. This implies that matter is not always compatible with F[R) gravity, not even in 
the Hofava-Lifshitz case. In our specific example po > implies 

X < U > U ^ < A < 6m^-3m 

m<0, 0<w<l, <(x>0/i<0 A> em^+"^-3w^-3^ 

^ ^ — 6m2— 3m 

= 0<m<i A< 6mM+m-3M 

2 2 2 ^™ 

ri ^ ]_ 12m fi+2ni — 15m/j— ^ \ ^ 67nfi+7n — 5fj, 

Y>0 U>0 i "^^2 6m^-3m A < ^ 

^ ' ^ '^0<ii;<10<m<i A< 6rnM+m-3^M-3A. ^^^"^ 

— 2 3m 
r\ ^ ^1 1^^^ l-\-w 6m^ii)u+12m^u+2m'^ — QmiLJu— 15m/j.— m+3iou+3u ^ \ ^ 6ma-\-m — 3iuii—3u 
U<W<i ^ <m< 6m2-3m <^< ^ 

Q m>i A < 12mV+2m^-15mp-m+3A. 

2 Grw^—Sm 
X<0, /i>0, ^0<W<li<m< -'-+"' A < 6'"'''"A'+12m^jj+2m''-9mM;/j-15m;j-m+3'M;j+3/J (141) 

< W < 1 m > i±2il'" A < 6'»A^+" -3^"3A' 

— 2i(j 3m 

ii; = m>i ^ 12mV+2m^-15mM-m+3M 

2 Gm'^— 3m 

Y>0 U<0 ^0<U'<l-<m< A < 6m^w^+12m^//+2m^-9mwM-15m^-m+3if;^+3M (142) 

^1 ^l+iu 6m'^i(j/^+12m^/j,+2m^ —9mwjj,— 15mfi—'m-{-3'Wjj,-\-3j-i , \ . 6mjj,-\-m — 3iufi—3j-i 
U<WSi > ^:^r 6m2-3m <A< 5;;^^ 

W = 0<m<i 12m''/x+2m''-15m/x-m+3M ^ ^ < 6mA.+m-3M 

2 bm^ —3m 3m 

<o <o<!'""'^ ""^^ """tm"'"'' n43) 

^ ' ^ ' In^ ^1 ^1 6rn^tuu+127n'^ u+2m'^ — 9miuu— 15m/j.— m+3iuu+3u . \ . 6mLL-\-m—3wLX^u ^ 
U<W<i U<?71<2 6m^-3m < < 3^;^ 

< w < 1 i < m < i±2^ A < s^m+^-^s^m^Sa. 

— 2 2u) 3m 

Another result which will be useful for our purposes is an exact solution for the theory F(R) = R + ^R? + xi?™. For 
TO > 2, this solution reads 

2 1 - 3A 

a = aot ' , 7 = — — , p, 



3(1 + w)' ^ 3(ii;-l)' 

2 4(3A-1) , 

P = Pot ' , Po = ^1^2 2 ■ 144 



This solution is obviously present only in the Hofava-Lifshitz version of this theory, as one can check directly. 
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One of the most important methods used to investigate power law solutions in higher order gravity is the recon- 
struction of a theory starting from a specific background. In the following we will adapt this technique to reconstruct 
the form of the function F{R) that admits flat FRW power law solutions 20j, |27H31|. 

Let us then consider a cosmological solution characterized by the Hubble parameter 

H^l, (145) 
and again assuming the energy density of a barotropic fluid 

P - Pot"^''^^+'"^ ■ (146) 

In this case the Ricci scalar is 

R ^ + + ^ (147) 

so that we can express the time < as a function of R. Substituting (|145p and (|146p into (|123p and (|126p and expressing 
t in terms of _R, one obtains 

AiR^F^^'> + A2R^F" + A^RF' + A^F + ewR^^'^+^^'' = , (148) 
BiR^F" + B2RF' + B3F + B4-R^('"+i)'^ = , (149) 

(150) 



with 



37 (7 - 37A + 67M - 2Ai) ' ^^^^^ 
_ 4(7-37A + 3m) 

37(7(3A-6m-1) + 2/x)' ^'^^^ 
2(37^1)(3A-3a.-1) 

^'-3(7-37A + 67M-2/.)' ^ ^ 

^4 = 3-i(^"+i)^Aj2u;po [7'(1 - 3A + 6^) - 2^ (154) 



and 



Bi = — ^^^^ , (155) 

37A — 67 /i — 7 + 2/i 

7(3A-1) 

7 — 37A + 67/i — 2^ 

S3 = 1, (157) 
^4 = «'po (-3-i('"+i)^) (7 (7 - 37A + 67/i - 2/.))-5('"+i)^ . (158) 

These equations admit the solution 

F{R) ^ CiR"- + C2^"+ + C3-R5(i+«')7 ^ (159) 

where 



and 



7 (3A -3^i-l) + 3^i± A/72 (-3A + 3Ai + 1)^ + 27/i (3A + 3^-1) + ^^ 

a± = ^- (160) 

4/i 



^ 3-i('"+i)'r^2^o [7 (-37A + 67M + 7 - 2/i)]^-^(-+^^^ 

^ 7[7(3A-l)(3(u; + l)7-l)-/x(3(u. + l)7-2)(3(w + 2)7-l)] ■ ^ ' 
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Note that the coefhcients a are real only for 

7^ (-3A + 3/i + 1)^ + 27^ (3A + - 1) + > , 

which is satisfied for 

7 > (162) 




Therefore also in the Hof ava-Lifshitz case, the only type of function F that is able to generate analytical power law 
solutions is a combination of powers of the Ricci scalar. The connection between the equations (jl45p and (|159p allows 
one to make some general considerations on the relation between the structure of the function F and the cosmology. 
If one plots the behavior of the exponents of (|159p as a function of 7 for various values of and A (see figure [1]), one 
finds that there is a correlation between the existence of a specific type of solutions and the sign of the modes of the 
reconstructed theory. For example, if one chooses only positive values for the exponents of (|159p in order to avoid 
instabilities, none of the permitted values of the parameters are able to generate a contracting solution. On the other 
hand, both a Friedmann expansion and a power law inflation regimes can be obtained if ^ < and A < ^"''^1^^^^^ 

or /i > and A > . The behavior of the exponents of (I159P is shown in Figure [T] for different values of the 

parameters. 

As a final remark, it is interesting to note that the case F[Ii) = C2R + CsR"^ corresponds to the solution (|136|) via 
the reconstruction method. This conclusion confirms the correctness of this approach and its utility in the search for 
exact solutions. 



C. Explicit model for the unification of inflation with dark energy 

It is interesting to try to formulate an explicit model for the unification of early-time inflation with late-time 
acceleration. One may propose such a model, which may unify the inflation and the late-time acceleration. First we 
consider the case 3A — 3/i — 1 = 0, which is specific in the Hofava-Lifshitz F{R) gravity. An example is 

F{R) = (1 + ci In K^R) (1 - C2K^R) . (164) 

Here ci and C2 are dimensionless positive constants. Then following (133), we find two de Sitter solutions 

1 

C2K 



R^ Rl = K^e ^ , R = Ri = — ^ . (165) 



If one chooses ci ~ 1/280, we find Rl ^ {iQ-^^ eY) , which may describe the accelerating expansion of the present 
universe. On the other hand, if C2 ^ C(l) ^ 0(100), R ^ Rj may express the inflation. 
In the general case 3A — 3/i — 1 7^ 0, we may consider the following form of F{R): 

/\/\ R R R R R R 

F i^Rj = R + f i^Rj f i^Rj = Ri tanh — + Rl tanh — + Ri tanh + Rl tanh -j- . (166) 

Here i?/, Rl, Ri, R2, and A are positive constants and we assume 

i?/>i?L>A i?/>i?i i?L>i?2. (167) 

Then, when R ^ Rj, we find 

f{R) - Ri , (168) 

such that f{R) plays the role of a large cosmological constant, which generates infiation. On the other hand, when 
R ~ Rl, f{R) becomes a small constant, 

f{R)^RL, (169) 
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(c)The exponents of l(T59)l {ov w = 0, fi = -2/5 and A = 1/2. 

FIG. 1: Plot of the curves representing the values of the exponents of the reconstructed F{R) theory corresponding to a 
background H = y/t, for u; = and some specific values of the parameters A and fi. From these plots one can infer, for 
example, that some backgrounds of this type can only be realized in F{R) theories with poles. 
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and the late time acceleration could be generated. Note that 

/(0) = 0, (170) 

therefore f{R) is not real cosmological constant. Hence, explicit construction of realistic models for the unification of 
the inflation with dark energy is possible. The remaining freedom in the choice of parameters gives the possibility to 
make the model quite satisfactory from the cosmological point of view. 



D. Finite-time future singularities 

The attempts of constructing accelerating cosmological models that include a dark component have revealed that 
such models often contain some unexpected phenomenology. One of the most striking features of dark energy cos- 
mologies is that, regardless of the way in which the dark component is introduced, they can become singular. By 
"become singular" we intend that there exists a specific time ts in which one of more key quantities of the model 
becomes divergent. Some of these singularities, like the so-called "Big Rip" 32|, are realized only far in the future (i.e. 
ts >> 1). However, as it was recently pointed out 33, S^], under special circumstances one can have quintessence-like 
cosmologies that present softer singularities at smaller finite time ("sudden singularities"). Although in pure GR 
cosmologies this pathological behavior can be cured by specifying, for example, the equation of state of the fluid, this 
is not the case in models that include dark fluids or in modifled gravity. In fact, it has been proved that these theories 
can admit up to four different types of singularities at flnite time [35j|. In the following, we will analyze the presence 
of these singularities in Hofava-Lifshitz F{R) gravity using some specific examples. 

Let us define the effective density and the effective pressure associated to the Hofava-Lifshitz F{R) gravity. We 
have 



1 



Pcff = 



Q^iHRF" - Q^iHF' + 6(3A - 3^ - l)H'^F' + F 
-2/iF(3)^2 „ 2(3A + /i - 1)HRF" - 2(3A - 3/i - 1)HF' - 6(3A - 3^^ - l)H'^F' - F 



(171) 



The different types of finite-time singularities can be classified by looking at the behavior of the quantities (|17ip plus 
the scale factor a, H and its derivatives. In particular, in they are classified as 

• Type I ("Big Rip") : For t ts, a ^> oo and (p, |p|) — > oo or p and \p\ arc finite; 

• Type II ("sudden") : For t tg, a ag, p Ps and \p\ — > oo; 

• Type HI : For t ^ tg, a ^ as, p ^ cxd and \p\ — ?■ oo; 

• Type IV : For t tg, a as {p, \p\) constant (or zero) and higher derivatives of H diverge. 

To classify the singularities in our case, let us consider the case of a vacuum spatially-fiat cosmology, and let us 
imagine that close to the time tg the Hubble parameter can be written as 



H^ho{ts-ty 



This means that the scale factor is 



a « tto exp 



hojtg - tf-i 
7-1 



(172) 



(173) 



The above expression tells us that there are two possible behaviors of the scale factor that depend on the value of 7: 
if 7 > 1, a will diverge as t approaches tg while if 7 < 1, a will converge. Therefore it is clear that the singularity of 
type I is realized when 7 > 1, the others for 7 < 1. 

The value of 7 also influences the form of the Ricci scalar. In general, one has 



R w hi (-9A + 18/i + 3) [tg - t)"27 _ GjhQfi{tg - ty'-''^^^ 



(174) 
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but, depending on the value of 7, the above expression can be reduced to 



^ ^ Mo (-9A + 18/i + 3) (t, - t)-^^ if 7 > 1 , ^^^^^ 



This property of the curvature also indicates that for t ^ one has 

f 1 7 > -1, 
I 1 7 < -1, 



(176) 



i.e. the curvature may become divergent or very small depending on the value of 7. This property will turn out to be 
very useful for simplifying the calculations. 

Finally, because of the nature of the structure of the effective energy density and pressure, theories with the different 
types of singularities are associated directly with the values of 7. Specifically: 

• Type I 7 > 1; 

• Type II =^ -1 < 7 < 0; 

• Type III ^ < 7 < 1; 

• Type IV ^ 7 < -1. 

Let us now consider some simple examples in the Hofava-Lifshitz F{R) gravity and their comparison with the standard 



1. The case F{R) = R + x-R" 

In the case 

F{R) ^R + xff" , (177) 

substituting (jl72l) and (I175p . one finds the necessary conditions for the presence of the singularities: 

Type I 7 > 1 

Type II TO < 0, -1< 7 < 
Type III < 7 < 1, TO 7^ 

Type IV 7 e Q - Z, 7 < -1 . 

These are compatible with the solutions found in [s^. It is important to stress that the conditions (I178P are only 
necessary for the existence of the singularity. The reason is that we have implicitly postulated that the solution (|172p 
satisfies (|123l) and (|126p at least in a specific time interval, which may not be the case due to the non-linearity of the 
theory. Then the only way to proceed is to find some exact solutions of the theory we are examining and see if the 
parameters have values for which the conditions (|178|) are satisfied. Unfortunately finding exact solutions in these 
theories can be problematic. However, close to the singularity one can solve this problem by using the fact that, as 
we have seen, the magnitude of the Ricci scalar also changes when t ^ tg. This means that we can approximate 
the function F with a simpler form that admits simple exact solutions. Then these solutions can be used to obtain 
necessary and sufficient conditions for the singularity to be realized. 
In our simple example it is clear that one has 




This means that we can approximate our F with i?" in all the cases of interest and that we can use the exact solutions 
found in Section FVlB I in order to understand the presence of singularities. In particular, one sees that the solution 



(|137p can present a singularity of Type I for 



A < 



a singularity of Type II for 



< TO < i 

1 < TO < 1 
1 < TO < 2 

TO = 2 
TO > 2 



6mA — 2m — 3A+1 ^ n ^ — 6mA+2m+3A — 1 
8mi'-16m^+16m-8 ^ P itn'^ -12m+8 ■ 

6mA — 2m — 3A+1 ^ ^ — 6mA+2m+3A— 1 

6m-6 <^ \i- <~ 4m^-12m+8 ' 

— 6mA+2m+3A — 1 ^ ^ 6mA — 2m — 3A+1 
4:m^ -12m+8 *^ M ^ 6m-6 ' 

/i< i(9A-3), 

^ 6mA-2m-3A+l ^ -6mA+2m+3A-l 

6m-6 4m^-12m+8 



TO > 2 



< TO < 5 

1 < TO, < 1 
1 < TO < 2 

TO = 2 
TO > 2 



6mA — 2m — 3A+1 



— 6mA+2m+3A — 1 ^ ^ 

4m2-12m+8 ^ A* ^ 8m3_i6m2+16m-8 ' 
-6mA+2m+3A-l ^ ,. ^ 6mA-2m-3A+l 
4m^-12m+8 ^ M ^ 6m-6 ' 



6mA-2m-3A+l 



M < 



6m — 6 <C *C 4777,^ 

/i> i(9A-3), 

-6mA+2m+3A-l 



6mA+2m+3A-l 



4m^-12m+8 



M > 



12m+8 ' 

6mA-2m-3A+l 
6m— 6 
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(179) 



(180) 



(181) 



TO < -1 < M < 4;;j7 

TO = — 1 /X > , 



6mA — 2m — 3A+1 
4 



-1 < TO < ^ < 
-1<TO<0 A^T^O 



6mA — 2m — 3A+1 
4m2 —4 



TO < —1 
TO = —1 
-1 < TO < /i < 



6mA — 2m — 3A+1 ^ , , ^ (\ 

T^^P^^ < M < 



Ai < 0, 

6mA — 2m — 3A+1 

M ^ 4m2-4 



(182) 
(183) 
(184) 



and a singularity of Type III for 



3 



TO < 

< TO < 

i<TO<l Ai<0 

1 < TO < 2 



^ 6mA — 2m — 3A+1 
A* 6m-6 ' 

-6mA+2m+3A — 1 ^ ^ (\ 
4m2-12m,4-S ^ A* ^ ^ : 



4m2-12m+8 

^ -6mA+2m+3A-l 
4m2-12m+8 



TO = 2 
TO > 2 

i < TO < 2 

2 

TO < 

< TO < i 

1 < TO < 1 
1 < TO < 2 

TO = 2 

TO > 2 



^ -6mA+2m+3A-l 
^ 'km'-' -12rn+8 

Ai > 0, 



< A* < 
TO 7^ 1 A« 7^ 



6mA+2m+3A-l 
4m^-12m+8 



Ai < M > fol-li^^^^ 
< Ai < 



4m^ 

6mA+2m+3A 



-6mA+2m+3A-l 



At < A* > 



-12m+8 
-1 



12m+8 A* > , 

6mA+2m+3A- 



4m^-12m+8 

Ai < 

6mA+2m+3A — 1 ^ ^ c\ 
4m2-12m+8 < A* < ^ ■ 



(185) 



(186) 



(187) 



This is very different from the standard case, where we have a singularity of Type I only for to > 2, and a Type III 
for i < TO < 1 and to < 0. 

For the solution (|136p we have a singularity of Type III when 



-3<TO<-| i(-2TO- 3) < u; < 1 
< TO < i (1 - x/T3) < w < 1 . 
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For the new solution (|138|) . which exists only in the Hofava-Lifshitz version of F{R) gravity, we have instead a 
singularity of Type I for 



and a Type III for 



f A < i i(3A - 1) < ^ < i(3A - 1) m > 1 
I A > i i(3A - 1) < ^ < i(3A - 1) m > 1 



A < i i(3A -1) < fi<0 m>l, 
A > i Q< fi< |(3A- 1) m > 1 . 



The results above show clearly that the presence of singularities is deeply altered in the Hofava-Lifshitz version of 
F{R) gravity. In particular, it seems that the additional parameters make it much easier to realize the singularities. 
The intervals we have presented above for the parameters can then be interpreted as constraints on this type of 
Hofava-Lifshitz F{R) theories of gravity. Thus, we have demonstrated that modified Hofava-Lifshitz gravity has the 
phantom-like or quitessence-like accelerating cosmologies, which might lead to singularities of type I, II, or HI. 



2. Eliminating the singularities 

Using conformal techniques, it was first argued in 36:] that in the standard F(R) gravity the singularities of the 
type we have found can be cured, when additional powers of the Ricci scalar are added to the Lagrangian. We can 
then verify if something similar happens also in the Hofava-Lifshitz case. Let us therefore consider the case 

F{R) =R + ^R^ + xi?" . (188) 

The action of this theory contains a correction of order R^ and, in the standard case, it is able to cure the singularities 
of theories of the type R + xi?™. 

If one derives the necessary conditions for the presence of a singularity, one finds: 

Type I 7 > 1, 

Type II never. 
Type HI < 7 < 1 m 7^ , 
Type IV jeQ-Z 7<-l. 

One can already see at this level that in this kind of theory singularities of Type II never occur. This can be interpreted 
as the fact that the correction R^ is able to compensate for these terms. 

As said before, the conditions above are only necessary and the only way to actually determine if a singularity is 
present is to consider an exact solution of the theory and see if the conditions above can be satisfied by that solution. 
Let us consider then the solution (jl44l) found in the previous section. Applying the conditions above one obtains 
that none of them is satisfied for this background. In other words, the addition of the R^ term compensates for the 
singularities one would find in a similar background of (jl77p . This supports the claim made in [36] that (jl88p is more 
regular than (jl77p and that in general the introduction of additional curvature invariants into the action can help 
in the curing of the singularities of an F{R) theory of gravity. On the other hand, it is also quite possible that the 



account of quantum gravity effects [37[ may also cure the future singularities. 

Summarizing, because of the additional parameters, the Hofava-Lifshitz version of F{R) gravity has a bigger space 
of de Sitter solutions compared to its standard counterpart. Using a simple theory and both a direct resolution of 
the cosmological equations and a reconstruction technique, we have also verified that this is true in the case of power 
law solutions. In general the presence of additional parameters also means that one has a bigger freedom in the 
choice of the features of these solutions, e.g. to see if it is possible to realize accelerated expansion. Therefore in 
these theories the realization of dark energy era (and inflation) is comparatively easier. This is interesting because 
it draws a direct connection between the Newtonian nature of quantum gravity and the observed behavior of the 
Universe. Unfortunately, the additional number of parameters also increases the probability that these cosmologies 
will become singular not only ed t ^ oo, but also at finite time. Using exact solutions, we have shown that there are 
many combinations of values of the parameters of the theory which are able to induce the appearance of singularities. 
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We have also verified, via a specific example, that adding an invariant of the type into the Lagrangian of a theory, 
one is able to obtain a theory whose solutions are much more regular. Consequently, also in the Hofava-Lifshitz case 
one can compensate such ill behavior of these models by the introduction of additional powers of the Ricci scalar into 
the action. 



VII. CORRECTIONS TO THE NEWTON LAW 



Let us now consider the possible corrections to the Newton law. For this purpose, we consider the infrared region 
where the higher derivative terms like (p4)) or (|39l) can be neglected and we find 



(g^^) ^ . (190) 

The term can be added from the beginning or this term might be induced at the infrared fixed point Then 
by the transformation (^5]) . by using ([55)1 and ([5n| . in (f^ . one gets 

d^xV¥^ |e^ _ ^g(3)^-'vf VvPV) - {A iv) F' {A (^))) - F {A (^)))| . (191) 

The usual Newton law can be generated through the exchange of the graviton. Furthermore by the exchange of the 
scalar field Lp, extra force might be generated. Now we consider the case that 

F{A) =A-Ks-^ + . (192) 

Here Aoff is an effective cosmological constant. Then 

P'iA) = ^ + ^,+---- (193) 

In the solar system or on the earth, the second term in (|193p is much smaller than unity, which corresponds to the 
first term. Hence, we find 



and therefore 



V{^) - Ae. + ^ Ae. + + . (195) 



Then since e*^ ~ 1, the effective mass rriip of if is given by 



3T1 + 2 



^2 = ^ ^^(2n+l)c^ ^ 3— (2n+l) 

o n+ I ^ „+i i)cTj+r 

In the "realistic" model, c is chosen to be c = ^^^^ ^ lO^'''^ cV. On the other hand, we find A = R ^ 

lO^'^^eV^ in the solar system and A ^ R ^ 10~^°eV^ on the earth. Thus, one finds ^ lO^^""^'' eV^ in the 
solar system and ~ lO'*^""'^^ eV^. Hence, if n could be large enough, the mass of ip would become large and the 
Compton length would become small, so that the correction to the Newton law would not be observed. 



VIII. DISCUSSION AND CONCLUSIONS 



We have proposed a first-order modified Hofava-Lifshitz-like gravity action and studied its Hamiltonian structure. 
As a large explicit class of such models we considered the modified Hofava-Lifshitz F{R) gravity that is more general 
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than the one introduced in ref. Q , which for the special choice of parameter fx = coincides with the degenerate model 
introduced in ref. [l5| . Its ultraviolet properties are discussed and it is demonstrated that such F{R) gravity may be 
renormalizable for the case 2; = 3 in a similar way as the original proposal for Hof ava-Lifshitz gravity. The Hamiltonian 
analysis of the proposed modified Hof ava-Lifshitz F{R) gravity shows that this theory is generally consistent with 
reasonable assumptions. The F{R) gravity action has also been analyzed in the fixed gauge form, where the presence 
of the extra scalar is particularly illustrative. The methods presented in the Hamiltonian analyses of sections IHII and 
IIVI can be used to study any action of the general form ([5]). 

The spatially-flat FRW cosmology of the modified Hofava-Lifshitz F{R) gravity is studied. It is shown that it 
coincides with the one of the earlier model but only in the spatially- flat FRW case. For specific choice of the 
parameters of the theory, its FRW equations of motion coincide with the ones of the traditional F{R) gravity. The 
presence of the multiple de Sitter solutions shows the principal possibility of the unification of the early-time inflation 
with the late-time acceleration in the Hofava-Lifshitz background, which proves that it can have rich cosmological 
applications. The power-law theories are investigated in detail. A number of analytical FRW solutions is found, 
including the ones with behavior relevant for the early/late cosmic acceleration. The quintessence/phantom- like 
cosmologies derived in our work may show all the four possible types of finite-time future singularities like in the case 
of standard dark energy. The conditions to cure such future singularities are discussed in analogy with the traditional 
F{R) gravity. It is also interesting that the correction to the Newton law in the F{R) gravity under discussion can 
be made unobservably small. Finally, a covariant proposal for F(R) gravity in Hofava-Lifshitz spirit has been made. 

Despite some successes in the formulation of modified Hofava-Lifshitz F{R) gravity which can be made renormal- 
izable and in its cosmological applications, a number of unsolved questions remain. What is the appropriate way 
to introduce matter in the theory? Is the theory itself fundamental (or at least, fully consistent) or does it descend 
from another more fundamental proposal? Can it comply with all the local tests in the Solar system as well as with 
cosmological bounds? What is the dynamical scenario for the restoration of the Lorentz invariance at late times? 
What are the cosmological and astrophysical consequences of the first-order modified Hofava-Lifshitz gravity when 
compared with those of the traditional modified gravity [39| . Moreover, the traditional questions about the properties 
of black holes in such a theory can be straightforwardly investigated. 

Nevertheless, even at the present stage some surprises can be expected from the theory. 

While the universe has likely undergone a perioid of inflation in its early moments, it is interesting to note that 
Hofava-Lifshitz gravity could produce cosmological perturbations that are almost scale- invariant even without inflation 

m. 



Hofava-Lifshitz gravity has also been considered in the presence of scalar fields |41l. |42| . In principle, it is possible 
to extend our Hofava-Lifshitz F{R) gravity by including its coupling with scalar fields. 

We would also like to mention a recent paper [43| , where a new class of Lorentz- invariance breaking non-relativistic 
string theories, inspired by the Hofava-Lifshitz gravity, has been presented and analyzed. 

Using the F(R) version of gravity one can propose even a more general formulation of string theory in the Hofava- 
Lifshitz background: for instance, rigid strings, membranes and p-branes, etc. On the other hand, it may suggest 
unusual solutions for the known cosmological problems. There also exists an attempt to explain the homogeneity of 
our universe in a model with varying speed of light A4\ . Having in mind that in the ultraviolet region the speed of the 
Hofava-Lifshitz graviton changes, one may speculate that the homogeneity of the universe may be described without 
the need for inflation. In any case, such a theory is both theoretically and cosmologically rich and it deserves further 
study. 
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Appendix A: Proposal for a covariant F{R) gravity 
In [s^l , a new type of covariant Hof ava-Lifshitz-like gravity has been proposed. The action is given, for z = 2 case 

(Al) 



s = J d*xV^ { ^ - " (^'^'^M- + /^™)' } 



and for z = 3 case 

S = J 
and for z = 2n + 2 case, 
5 = 



R 

2^ 



a {T^^R^, + I3TR) (T^'^V^V, + -fTVPVp) {Tt^'R^, + pTR) 



(A2) 



(A3) 



Here T^i^ is the energy-momentum tensor of the perfect fluid with constant equation of state parameter w and the 
parameters /3 and 7 are given by 



w — 1 
2(3iy- 1) 



7 



3w- 1 



(A4) 



Note that this fluid is not the standard matter fluid. It may have a stringy origin or it may be a kind of gravitational 
fluid. Hence, we may consider the following covariant Hofava-Lifshitz-like gravity: 



S 



^ J d^x^FiR^o.) , 

R - 2aK^ (T'^'^R^^ + PTRf , z = 2 , 

R - 2K^a {Tf^'^R^, + pTR) (T^'^V^.V, + jT\/P\/p) (T^'^i?^, + PTR) z = 3 , 

R - 2K^a {Tf^^Rp,^ + PTR) [Tt^'^V^.V^ + jTVPVp) {T^^R^,^ + I3TR) z = 2n + 2. 



(A5) 



The theories given b y th e action (|A5|) could be renormalizable when z > 3. This may be demonstrated using the same 
arguments as in ref. [38| . This is a large class of covariant F{R) gravities whose ultraviolet properties are more similar 
to the ones of Hofava-Lifshitz-like gravity. However, in many respects their cosmology is similar to the cosmology of 
traditional modified theories of gravity [1| . 
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